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Abstract. Direct photon emission in heavy-ion collisions is calculated within the relativistic microscopic 
transport model UrQMD. We compare the results from the pure transport calculation to a hybrid- model 
calculation, where the high-density part of the evolution is replaced by an ideal 3-dimensional fluiddynamic 
calculation. The effects of viscosity, present in the transport model but neglected in ideal fluid-dynamics, 
are examined. We study the contribution of different production channels and non-thermal collisions to 
the spectrum of direct photons. Detailed comparison to the measurements by the WA 98-collaboration are 
undertaken. 

PACS. XX.XX.XX No PACS code given 

1 Introduction hint to the creation of a strongly coupled QGP (sQGP) [3, 

4,5,6]. Possible evidence for the creation of this new state 

The creation and study of high-density and -temperature of matter has also been put forward by collaborations at 

nuclear matter is the major goal of heavy-ion experiments, the Super Proton Synchrotron (CERN-SPS), as for in- 

If the energy density in the reaction is high enough, a state stance the step in the mean transverse mass excitation 

of quasi-free partonic degrees of freedom, a Quark-Gluon- function of protons, kaons and pions and the enhanced 

Plasma (QGP) [1,2], may be formed. Recent experimental if^/7r+-ratio [7]. 

observations at the Relativistic Heavy Ion Collider (BNL- Unfortunately, it is not yet possible to describe the 

RHIC) like e.g. strong jet quenching and high elliptic flow time evolution of the produced matter from first principle 
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Quantum Chromo dynamics (QCD). Neither can the hot is sufficiently dense, so in the late stages of a heavy-ion 

part of a collision be observed directly. Therefore well- collision fluid-dynamics looses applicability. An advantage 

developed dynamical models to describe the space-time is that in the dense stages, hydrodynamics can propagate 

evolution of nuclear interactions are needed. Among those any kind of matter, and also allows for transitions be- 

approaches is relativistic transport theory [8,9,10,11,12, tween two types of matter, e.g. QGP and hadron gas, if 

13]. In this kind of microscopic description, the hadronic an appropiate Equation of State (EoS) is provided. Fluid- 

and/or the partonic stage of the collision can be described dynamics can therefore be used to study hadronic and 

under certain approximations. E.g. most transport mod- partonic matter in one common model. 

els cannot describe collisions with more than two incoming T^e restrictions of this kind of model can be loosened 

particles restricting the applicability to low particle densi- as well. By introducing viscosity and heat conductivity, 

ties, where multi-particle interactions are less important, perfect thermal equiUbrium does not have to be present 

Some attemps to include multi-particle interactions do ex- ^ny point. However, even with second order corrections 

ist [14, 11, 15, 16, 17], but this field of study is still rather ^i^e matter has to be close to equilibrium. 

new. The coupling of a partonic phase with a hadronic ^ , , • , • , • , 

Input to solve the hydrodynamic dirterential equations 

phase imposes one more problem on transport models. , , . . . ... , ^ ^ , .i ,. 

are the boundaries, i.e. the initial state (the distributions 

The microscopic details of that transition are not known. r , j_ , i j_ j_ j_- i 

01 all relevant densities and currents at the time the evolu- 

There are several attempts to address the problem of mi- . , \ c c^, , ■^■ ,^ 

tion starts), the Equation of State providing the pressure 

croscopic hadronization [12,18,19,20,21,22]. One caveat r r- .i i i i •,• .i . i 

as tunction oi the energy and baryon densities, that dc- 

in transport approach is that all microscopic scatterings ., , , . r.^ .i . • • i ^ ^ 

scribes the behaviour oi the matter that is considered, and 

are explicitly in the model and therefore the cross-sections j. ^ , r 

the freeze-out hypersurtace. 

for all processes must be known or extrapolated. However, 
for most processes no experimental data are available, and 

Out of the many possible observables, electromagnetic 

a large fraction of cross-sections have to be calculated or 

probes have the advantage of leaving undisturbed: once 

parametrized by additional models. 

they are created, they will escape the reaction zone, due 
Relativistic, non-viscous, perfect fluid-dynamics is a to negligible rescattering cross-sections. Besides dileptons, 
different approach to explore the space-time evolution of direct photon emission is therefore of greatest interest to 
a heavy-ion collision. It constitutes a macroscopic descrip- gain insight into the early, hot and therefore possibly par- 
tion of the matter that is created, assuming that at every tonic stages of a reaction. Direct photons are distinguished 
time and in every place it is in perfect local thermal equi- from the mass of photons as those coming from collisions 
librium. This assumption can only be true if the matter and not decays. 
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Unfortunately, most photons in heavy-ion coUisions eraged over in order to apply thermal rates to the coarse- 
come from hadronic decays in the very late stages, mostly grained distributions [23]. 

1^0 _^ These decay-photons impose a serious chal- The application of microscopic cross-sections can only 

lenge to the experimentalists when they try to obtain the be undertaken in a model where all microscopic collisions 

spectra of direct photons. Up to now, several experiments are known. That limits the field of use to transport models, 

have tried to obtain the spectra of direct photons: He- There have been several calculations for photon spectra 

lios, WA 80 and CERES (all at CERN-SPS) could pub- from transport models [24,25] 

lish upper limits, while WA 98 (CERN-SPS) and PHENIX In this paper we first introduce the dynamic mod- 
(BNL-RHIC) have published data points for direct pho- els used for the calculations. Then, we explain the cross- 
tons, sections and thermal rates used for photon emission as well 



Approaches to the theoretical description of direct pho- 
ton spectra include calculations with perturbative Quan- 
tum-Chromo-Dynamics (pQCD), the application of ther- 
mal rates and microscopic cross-sections, all of which have 
their area of applicability. 



as the mechanisms for doing so. Following that, the results 
of our calculations are presented and compared to results 
from the WA 98-collaboration. Finally, after summary and 
conclusions, an outlook to further planned studies is given. 

2 UrQMD hybrid model 



Calculations based on pQCD describe the photon data 
in proton-proton collisions very well and, if scaled by the 
number of binary collisions, also those in heavy-ion reac- 
tions. The range of applicability is limited to high trans- 
verse momentum p]^^ 1 GeV. 



UrQMD (Ultra-relativistic Quantum Molecular Dynam- 
ics) is a microscopic transport model. It includes all ha- 
drons and resonances up to masses m = 2.2 GeV and at 
high energies can excite and fragment strings. The cross- 
sections are either calculated via detailed balance or para- 
Thermal rates can only be applied if the assumption metrized by the additive quark model (AQM), if no ex- 
of local thermal equilibrium is fulfilled. Scattering rates perimental values are available. At high parton momen- 
can then be calculated by folding the particle distribu- tum transfers, PYTHIA is employed for pQCD scatter- 
tion functions of the participating particle species with ings. UrQMD has been used by Dumitru et al to study 
the respective cross-sections. This framework can be ap- direct photon emission earlier [24] . 

plied to either static models, simplified hydrodynamics- In this work, we combine and compare the two models 
inspired models such as the blast wave model and to full mentioned above to describe the space-time evolution of 
fluid-dynamic calculations. The space-time evolution of a a heavy-ion reaction: for the unequilibrated initial state 
reaction as predicted by microscopic theories can be av- and the low-density final state, the microscopic transport 
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model UrQMD is used, whereas the high-density part of 3 PhotOtl emission from hadronic SOUrces 

the reaction is modelled using ideal 3+1-dimensional fluid- 

In both models, hybrid and pure transport, photon emis- 

dynamics. For details of this model see [26]. During the 

sion is calculated pcrturbativcly. That means, the evolu- 

evolution, two mappings have to be performed: After the 

tion of the underlying event is not changed by the emission 

incoming nuclei have passed each other, the baryon-number- 



, energy- and momentum-densities are smoothed and put 
into the hydro calculation. After the local rest frame en- 
ergy density has dropped below a threshold value of ecrit = 
730 MeV/fm^ (w 5eo) in every point, particles are created 
from the densities by means of the Cooper-Prye formula 
and propagation is continued in UrQMD. The effect of 
changing those parameters is studied in [26] and found to 
be small. 



of photons. This is justified as the emission probability for 
photons is extremely small. 

3.1 Photons from UrQMD 

In the transport part of our calculations for each scatter- 
ing the cross-section for photon production is calculated. 
Cross-sections are taken from Kapusta et al [27] and Xiong 
et al [28] . In the former work, the photon self-energy from 
a Lagrangian involving the pion, rho and photon-fields has 
been taken as basis for the calculations: 

For the present studies, the EoS used in the hydro- 
dynamic part resembles a free hadron gas with the same £ = |£)^#|^ — to^|#|^ — -p^ivP^" + -"m^pPixP^ — -FfivF^" 
degrees of freedom that the transport part (i.e., UrQMD) (1) 
has. Since strings are expected to be created in scatter- (for details the reader is referred to [27]). The coupling 
ings with very high center-of-mass energies, and since hy- of the p- to the 7r-meson is characterized by the couphng 
drodynamics only describes soft scatterings, they are not constant Qp, which is inferred from the decay rate p ^ tttt 
included in the EoS. with qI/Att = I2rujl/{mppl). Here, po and ujq are the mo- 

mentum and energy of a decay- tt in the rest frame of the 



By this choice, one can study the impact of different 
underlying dynamics of the reaction on the spectra of pho- 
tons. This allows to disentangle which changes are due to 
different physical assumptions. In further work, we plan 
to explore to the EoS, i.e, changes to the assumptions of 



p and r is the decay-width. Kapusta et al argue that F is 
experimentally measured and therefore effectively includes 
certain higher-order effects such as vertex corrections. 

Prom this Lagrangian, they calculate cross-sections for 
the following processes: 



the underlying medium (chirally restored or deconfined ~^ ' 

matter). tt^tt^ — > 7P^, 
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7r^p° — > 77r^, sion instead of only one (complete) photon every 10,000 

^±pT _^ j^o^ collisions is created. This allows us to have a much bet- 

n + + ter statistics with less hadronic events. The small cross- 

n p 77r , 

sections for the processes we consider justiiy this pertur- 

n TT^ — > 'yy. 

bative ansatz. 

The last of these is suppressed by an additional factor 
of a with respect to the others and is therefore not ex- 

3.2 Photons from hydrodynamics 

pected to contribute to the photon spectra significantly. 

Kapusta et al also quote cross-sections involving the rj- For the hydro-evolution we also produce fractional pho- 
meson (tt^tt^ — > 77/ and tt^t? — > 777^), but they are omit- tons at every cell of the hydrodynamic calculation. We 
ted in our calculations. use the parametrizations by Turbide, Rapp and Gale [29]. 

From the paper by Xiong et al, we deduce the cross- To obtain these, they start with an effective non-linear a- 
sections for the processes model Lagrangian with the vector and axial vector fields 

implemented as massive gauge fields of the chiral U{3)l'x 

Trp ^ ai ^ 77r, 

U{3)r symmetry. They calculate a different set of rates 
averaged over all possible charge combinations. This chan- ^Yian cross-sections calculated by Kapusta et al, namely: 
nel is not included in Kapusta et al. 

TTTT 7p, 

In order to obtain photon spectra, all scatterings that 
happen during the transport phase are examined. For ev- ~^ ^''^^ 

ery scattering that can produce photons, the correspond- ^-^* ~^ 7-^) 

ing fraction of a photon ttK ^K*, 

j^t ^ ^ ^2) ^ 'y^' 

K*K 77r. 



is produced. In this formula, (Ttot is the sum of the total 

hadronic cross-section for a collision with these ingoing The rate for -Kp — > 77r contains the process np ^ ai ^ 'yn. 

particles (known from UrQMD) and the electromagnetic The rates involving K and K* are included in our calcu- 

cross-section a^m as calculated by the formulae from the lation albeit they are negligible compared to the others, 

abovementioned papers. Turbide et al also calculate the emission rate from p-decay 

This number of photons is then distributed in the {p — > inr'y), but this is not taken into account for the 

solid angle by the angular distributions as given by da/dt, present work, because we want to restrict ourselves to col- 

so that a whole distribution of fractional photons per colli- lisional photons. 
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For every cell, we generate one fractional photon, its 
fraction being the integral over the invariant rate: 



(3) 



where AV and At are cell volume and time step, respec- 
tively, and EdR/(Pp is the rate for photon emission as 
given by Turbide et al. 

To get photon spectra with the right lab-frame distri- 
bution out of the (spherically symmetric) local-rest-frame 
distribution, a lorentz-transformation has to be applied to 
the generated photons. 

We choose the invariant generalization of the energy 
p^u'^ according to the distribution functions implied by 
the rates. Then we choose a which, together with the 
cell's flow velocity m'', results in the desired p^u^ . 

One fractional photon is then created in the direction 
of this p^^ with the fraction N}^. 

4 Results 

All calculations are done for Pb+Pb-coUisions at incident 
beam energy of ii^Lab — 158 AGeV. The sample of col- 
lisions has impact parameters of 6 < 4.5 fm and only 
midrapidity-photons (|yc.m.| < 0.5) are included in the 
analysis. These settings are equivalent to the WA 98 trig- 
ger conditions and detector setups for their "central" data 
set [30,31]. 

4.1 Emission times 

Due to formation time effects and the Heisenberg princi- 
ple, particles cannot be produced earlier than tprod ~ E^^ . 
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Fig. 1. (Color Online.) Average emission time of direct pho- 
tons vs. transverse momentum. The curve for tttt 77 is omit- 
ted from the plot; the emission times for that rate are much 
higher than the range of the plot. 

For photons at mid-rapidity, where E = p±, this gives a 
lower bound for particle production times in different p_L- 
bins. 

Later emission times, though, are well possible. If a 
particular scattering occurs preferably at late times, the 
average emission time of photons may be shifted to later 
times. 

In this analysis only scatterings between mesons are 
considered as photon sources. Mesons have to be produced 
in a heavy ion collision, so it can be assumed that meson- 
meson scatterings occur only after some time. 

Figure 1 shows the average emission times of all pho- 
tons as well as of those photons coming from -np- and 
TTTT-scatterings. The complete average is obviously domi- 
nated by the processes with tt and p in the initial state, 
as is consistent with the results in section 4.2. 
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It is eye-catching that the average emission times of the In UrQMD, the leading particles from a string have a 

photons deviate a lot from the hyperbolic behaviour men- reduced cross-section during their formation time. The ef- 

tioned earlier. Especially in the intermediate p_L-region of fects of neglecting photons coming from collisions of string 

1 — 2.5 GeV it becomes apparent that photons are pro- ends is shown in Figure 2. The difference becomes impor- 

duced from decay products that have not been present in tant at high p±, where photons from initial hard pQCD- 

the collision before 5 or 6 fm. TTTr-scatterings contribute scatterings play an important role, 
at the same times, but at lower momenta {p± ~ 0.5 — 

1.5 GeV). Pqj. comparison, we show the pQCD-spectrum in Fig- 



4.2 Photon spectra 



ure 2 (taken from [33]) where an intrinsic {k'j_) « 0.9 GeV^ 
was used. One clearly observes that pQCD photons are 
a dominant source of photons at high p±, especially if 
the (slightly artificial) contribution from string ends is re- 
moved from the spectrum. 



Figure 2 (left) shows the direct photon spectra from a pure 
transport calculation. The most dominant contribution to 
the spectra is the channel irp jtt, which contributes 
about on order of magnitude more than the sum over all 

other channels between = 0.5 GeV and p± = 2.5 GeV. ^^'^ » comparison between hybrid- and transport-cal- 
Emission from all channels is very thermalized at low culation we restrict the analysis to those channels that 
transverse momenta, and the nn ^ 7p- and Trp ^ jn- ^rc present in both models. Figure 3 shows spectra for 
channels contribute significantly at high p_L, where pre- ^ TP plus np jtt for pure transport and the corn- 
equilibrium scatterings play an important role. P^^te hybrid model. The different stages to the hybrid 
As expected, the doubly-electromagnetic channel tttt ^ "^o^el (pre-equilibrium, hydro, post-freeze-out) are seper- 
77 is very much suppressed over the whole p_L-range. ately plotted. 

The relative contributions of the different channel are 

very similar in the case of a hybrid-model calculation. Fig- The yield at low p± < 1 GeV is dominated by radiation 

ure 2 (right) shows the spectra obtained in that case. Here, from the hydro phase. Above that value, most photons 

the dominance of the 7rp-channcl is as pronounced as in come from the post-freeze-out and pre-equilibrium stages, 

the pure transport calculation. Overall, contributions from It is interesting to see that the photon yield in the range 

that channel and the TTTr-channel look rather similar to the 1 GeV < p± < 2.5 GeV is dominated by photons coming 

UrQMD-case above a momentum of about p_L « 1 GeV, from post-freeze-out scatterings. This is consistent with 

because in this regime the non-equilibrium phases con- the observation of late emission times in that p_L-range as 

tribute more than the high-density hydro phase. shown in Figure 1. 
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Fig. 2. (Color Online.) Transverse momentum spectra of direct photons from both models. Left: pure UrQMD, right: Hybrid. 
The different charge combinations of the channels have been combined for better clarity. Experimental data taken from [30] 
("central" data set, table IV), arrows indicate data points where the error band extends to negative yields. In the left figure, 
a pQCD-spectrum from the initial proton-proton-collisions [32,33] is shown, as well as the spectrum from pure UrQMD when 
neglecting the contribution from string ends (see text). 

used and compared two different models, both of which 
give similar results. UrQMD, as a purely hadronic trans- 
port model, underpredicts the data from WA 98 below 
p± — 3.5 GeV. The hybrid calculation, where the high- 
density phase of the evolution is described by hydrody- 
namics, also underpredicts the data in that range. 

In both cases, processes with n and p in the initial 
state dominate the intermediate and low momentum re- 
gions and processes with two tt become important at high 
momenta, where early scatterings dominate. 

Emission times for direct photons at intermediate trans- 
verse momenta are found to be very late, and so this p±- 
range is in the hybrid model dominated by photons from 
the late post-freeze-out stage. 
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Fig. 3. (Color Online.) Transverse momentum spectra of direct 
photons compared between both approaches. 

5 Summary 

In this work, we have shown the direct photon spectra for 
central Pb-|-Pb-collisions at top SPS energies. We have 
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At very low p±, the hydro phase contributes much the r/- meson as well as baryonic processes should be among 
more to the photon spectra than the transport phase does, those that are added next. 

Here, the only significant difference between both models Also, the existing channels can be improved: Until now, 
is seen. the channels that involve a p-meson as outgoing particle 

always assume that the p-meson is produced at its pole 
mass, which clearly does not have to be the case. 

6 Conclusions 

On the part of the hydro evolution, it is of course of 

The hadronic models used in the current work cannot re- ^'^^^ '""^^'^^^ *° ^^^^^ °f d^"^^^* P^°*°^« 

jj-uj-j-ui 4. ju iirA no thermal partonic processes, i.e. the inclusion of an EoS 

produce the direct photon spectra as measured by WA 98. ^ ^ ' 

rri, r 4.U- • 1 i u u ii, i ii, containingaphasetransitiontoaQGPisnecessary.com- 

ihe reason for this remains unclear; it could be that the ^ t- j 

^ r 1, 1 • Tj. • 1 parisons to data from RHIC and predictions for the new 

implemention or more channels is necessary, it is also pos- ^ ^ 

sible that using an EoS with phase-transition to a QGP FAIR-facility will be accessible then, 
will yield significantly higher spectra. The addition of pho- 
tons from hard pQCD-scatterings will contribute signifi- 8 Acknowledgements 
cantly, though not sufficiently, to the final spectrum. Most 



probably, each of these points will contribute to the miss- 
ing photons. 

7 Outlook 



This work has been supported by the Frankfurt Center for 
Scientific Computing (CSC), the GSI and the BMBF. The 
authors thank Hannah Petersen for providing the hybrid- 
and Dirk Rischke for the hydro-code. B. Bauchle grate- 
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The work presented here lays the foundation for many fur- 
ther investigations. First, photons from initial hard pQCD- 
scatterings have to be added. Other collisions with high 
^/s have to be treated with pQCD instead of hadronic 
cross-sections as well. The channels that are implemented 
in both stages of the model (transport and hydrodynam- 
ics) have to be aligned, so that a comparison does not ^ j ^^ Harris and B. Muller, Ann. Rev. Nucl. Part. Sci. 46 
have to be confined to the smallest common denominator. (1996) 71 [arXiv:hep-ph/9602235]. 

Also, the list of channels in both cases has to be extended 2. S. A. Bass, M. Gyulassy, H. Stocckor and W. Greiner, J. 
beyond what is implemented so far. Scatterings involving Phys. G 25 (1999) Rl [arXiv:hep-ph/9810281]. 
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